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Abstract 

The most recent development of the quark-meson couphng (QMC) model, in which the effect 
of the mean scalar field in-medium on the hyperfine interaction is also included self-consistently, 
is used to compute the properties of hypernuclei. The calculations for A and H hypernuclei are of 
comparable quality to earlier QMC results without the additional parameter needed there. Even 
more significantly, the additional repulsion associated with the increased hyperfine interaction in- 
medium completely changes the predictions for T, hypernuclei. Whereas in the earlier work they 
were bound by an amount similar to A hypernuclei, here they are unbound, in qualitative agreement 
with the experimental absence of such states. The equivalent non-relativistic potential felt by the 
S is repulsive inside the nuclear interior and weakly attractive in the nuclear surface, as suggested 
by the analysis of S-atoms. 
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I. INTRODUCTION 



The study of A hypernuclei has a long and impressive history, with the shell structure 



mapped out across the periodic table 



! 2, ^, 4, 5|. It is known that the single particle spin- 
orbit force is very small 6| and systematic studies of the energy levels of light A hypernuclei 
have enabled the extraction of considerable detail concerning the effective A-N interaction. 
Current studies of electroproduction at JLab will provide important new information in this 
area j^, 0] • 

When one turns to E and S hypernuclei, the situation is quite different. The special case of 
^He aside, there is no experimental evidence for any S hypernuclei s], [o], Q , despite extensive 
searches. Indeed, it seems likely that the S-nucleus interaction is somewhat repulsive and 
that there are no bound S hypernuclei beyond A=4. In the case of the H, the experimental 
situation is very challenging, with just a handful of observations of doubly strange nuclei. 
While we eagerly await studies of S hypernuclei with new facilities at J-PARC and GSI- 
FAIR, it seems likely that these hypernuclei may exist UJ, ll2l] and it would be helpful to 
have a range of predictions from various theoretical models. 

There is currently considerable interest in the equation of state of dense nuclear matter, 
especially in connection with the calculation of neutron star properties [isl, Q] • The density 
at which hyperons appear and particularly their order, plays a crucial role in the equation of 
state (EoS) - for example, the inclusion of hyperons in the QMC model seems essential if the 
direct URCA process is to play a role in cooling for observed neutron stars [l^. For matter 
in /^-equilibrium, negatively charged hyperons are favoured a priori and in many models S~ 
hyperons enter soon after the A. Other models suggest that the S~ should appear after 
the A (followed by the S*^), with the S playing no role. Whichever scenario is ultimately 
correct, hyperon physics will play a critical role at 2 - 3 po and a relativistic treatment is 



essential 



17|. 



In this work we apply the latest development of the quark-meson coupling (QMC) 
model to calculate the properties of hypernuclei. The major improvement in the model 
that we use here is the inclusion of the effect of the medium on the hyperfine interaction. 
This has the effect of increasing the splitting between the A and S masses as the density 
rises. This is the prime reason why we find that S hypernuclei are unbound. On the other 
hand, the model produces quite good binding energies for A hypernuclei and predicts modest 
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binding of S hypernuclei. 



II. RECENT IMPROVEMENT IN THE QMC MODEL 



The QMC model was created to provide insight into the structure of nuclear matter, 



starting at the quark level 



18 



19 



20l |. Nucleon internal structure was modeled using the 



MIT bag, while the binding was described by the self-consistent coupling of the confined 
quarks to a and uj meson fields generated by the confined quarks in the other "nucleons" in 
the medium. While the use of effective scalar and vector fields to carry the forces is very 
similar to QHD, the explicit treatment of the N internal structure represents an important 
departure from that classical approach. The self-consistent response of the bound quarks 
to the mean a field leads to a novel saturation mechanism for nuclear matter, with the 
enhancement of the lower components of the valence Dirac wave functions, as the density 
increases, reducing the effective a-N coupling. This effect is naturally represented by a single 
new parameter, d, the scalar polarizability, with the nucleon effective mass taking the form: 



M 



N 



(1) 



The QMC model has been used to study the binding of u, rj, rj' and D nuclei 2l|, l22|, |23|], as 
well as the effect of the medium on and J/\l/ production 
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. It leads to a very natural 
A recent extension of the 



explanation of the small spin orbit force in A hypernuclei 125 
same physical ideas to a confining version of the NJL model |27j has led to a quantitative 
description of the EMC effect across the periodic table, with a remarkable prediction of a 
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29|. 



sizeable enhancement of the nuclear modification of spin structure functions 

For the present purpose, the most significant recent development is the inclusion of the 
self-consistent effect of the mean scalar field on the familiar one-gluon-exchange hyperfine 
interaction that in free space leads to the N-A and S-A mass splitting 15|. The original 
QMC model treated this term as a small correction and ignored its medium modification. 
However, the term "hyperfine" is misleading, as the N-A splitting is 30% of the nucleon mass. 
One-gluon-exchange in the bag model is essentially determined by the magnetic moments 



of the confined quarks 



30| and, in first approximation, the hyperfine splitting is 



3 AT 



(2) 
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Tils 


A 


S 




S* 




n 


Fs = l 


0.341 


1.135 


1.176 


1.355 


1.416 


1.599 


1.784 


Fs = 0.726 


0.297 


1.107 


1.189 


1.325 


1.368 


1.507 


1.654 


Exp. 




1.116 


1.195 


1.315 


1.385 


1.533 


1.672 



TABLE I: Octet and decuplet masses in GeV. The parameter Fs and the strange mass nis are 
fitted to the A, S,H masses. The last three columns are predictions for the decuplet masses. 

where C is a constant, essentially independent of the quark masses, and /io,/is are, respec- 
tively, the magnetic moments of the confined light and strange quarks. From the experi- 
mental value in free space (5e — 6a)/ (5a — 6n) = 80/300 we get /is//io ~ 0.6. This decrease 
of the magnetic moment with the increase of the quark mass is a relativistic effect similar 
in nature to the saturation mechanism of the QMC model. In the medium the light quark 
interacts with the nuclear a field and this amounts to decreasing its mass. This in turn 
reduces its eigenenergy (which is roughly the equivalent of the constituent quark mass) and 
as a result the magnetic moment of the light quarks in medium changes to (1 + e)/io with 
e > 0. An explicit calculation predicts e ~ 7% at nuclear matter density. Assuming that 
the strange quark is not coupled to the nuclear field one thus finds that the relative change 
of 5s — S\, because of the a field, is 

(^e-^a)^ (l + e)[/io(l + e) -/i,] 



3.5e, 



(3) 



where we have used /is/ yUo ~ 0.6. This rough calculation thus implies that the S will be less 
bound than the A by 3.5e(5s — 6a) ~ 20MeV, which is consistent with the exact calculation 
described below. Our motivation for applying the model to finite hypernuclei is that this 
effect has the right sign to help solve the major problem found in the original QMC model, 
namely that S hypernuclei were not much less bound than A hypernuclei. 

In view of the physical importance of the hyperfine interaction, as we just explained, we 
decided that it would be worthwhile to improve on the MIT bag model, which traditionally 
yields only 50 MeV of the S — A mass difference, instead of the experimental value of 80 
MeV. The change which we introduce is motivated by the work of Barnes 3l| , who explained 
the relatively large value of the effective strong coupling constant, «<j, required in the naive 
bag model. That model ignored the enhancement of the relative q-q wave function at short 
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distance caused by the attractive color Coulomb force. The short distance enhancement 
of the relative wave function is expected to be less effective for heavier quarks. As a phe- 
nomenological way of implementing these effects, we multiply the hyperfine interaction in 
the ususal MIT bag model by a factor F", with n the number of strange quarks participat- 
ing. We fix the bag model parameters in two steps. First, the masses of the nucleon and A, 
together with their stability conditions, are used to determine the vacuum pressure, -B, zero 
point energy, Zq, and the strong coupling, a^, for a given value of the nucleon radius, which 
is set to 0.8fm in practice. This step is independent of the quark mass and F^. The nucleon 
and A masses are exactly reproduced with B = 0.5541 fm~^, Zq = 2.6422 and a.^ = 0.4477. 
Then the strange mass, m^, and the parameter Fg are chosen to give the best fit for the 
A, S and S masses. As shown in Table [U where we also display the masses corresponding 
to Fs = 1, we obtain an excellent fit with Fg = 0.726 and = 0.297 GeV. In particular, 
the A — S splitting is well reproduced. For completeness, we show in the last 3 columns 
of Table [T] the good agreement obtained for the predicted masses for the rest of the baryon 
decuplet (S*, S* and Q). 



III. APPLICATION TO THE ENERGY LEVELS OF HYPERNUCLEI 



In order to calculate the properties of finite hypernuclei, we construct a simple shell 
model, with the nucleon core calculated in a combination of self-consistent mean scalar 
and vector mean fields. The determination of the scalar coupling to the octet baryons 
requires a sophisticated self-consistency calculation. However, within the Born-Oppenheimer 
approximation, which requires that the internal structure of the nucleon has time to adjust 
to the local environment (an approximation estimated in Ref. 19| to be good at the level of 
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3% or better in finite nuclei), the result can be parametrized in a very practical form 151 ]: 



M. 



N 



MM 



+ 
Ma 
+ 

Ms 
+ 
M= 
+ 



gacr 

0.002143 + 0.10562/?^''"" - 0.01791 ( R{ 



•-TV 



0.6672 + 0M638Rj^^' - 0.0022 ( 



(4) 



0.00146 + 0.0691i?j^''^' - 0.00862 ( 



0.6653 - 0.08244i?{''^ + 0.00193 ( i?^''^ 



2" 



(5) 



0.00064 + 0.07869i? 



free 
N 



0.0179 ( R^^^' 



0.3331 + 0.00985i?J^ 



■AT 



-0.00032 + 0.0388i?^''^' 



0.00287 ( Rf;^""' 



0.0054 ( i?^'"" 



2 



2 



(6) 



(7) 



We take the bag radius, Rj^'^'^, to be 0.8fm but note that (c.f. Fig. 1 of Ref. [15|) the results 



are quite insensitive to this parameter. Note that the coefficients in Eqs. (jlj) to ([7]) differ 
slightly from those in Ref. 15|] because of the improvements in the MIT bag model explained 
in Sect. II. The only significant change is for the S hyperon and this is a result of ensuring the 
correct hyperfine splitting from the A. It is worth recalling that one underlying hypothesis of 
the model is that the (a, u, p) mesons do not couple to the strange quark, which is justified 
by the OZI rule since the interaction cannot proceed by quark exchange. Moreover, this 
hypothesis directly leads to the absence of spin-orbit splitting for A hypernuclei, which 
is consistent with experimental data. In this framework, the cj-baryon coupling goes like 
the number of non-strange quarks and thus follows the constraints of SU(6) symmetry. 
We note that, for example, QCD sum-rule investigations |l^ have suggested that the uj- 
baryon couplings may differ from the SU(6) values. On the other hand, the success of the 
QMC model in leading to very realistic effective NN forces gives us some confidence, a 
posteriori, for the hypothesis. 

To calculate the hyperon levels, we use a relativistic shell model. In principle it would be 
preferable to generate the shell model core, consisting only of nucleons, using the Hartree- 
Fock approximation. However, in practice it is much easier to use Hartree approximation, 
which should produce very similar results, because the coupling constants are adjusted in 
each case to reproduce the properties of nuclear matter. Thus, we use the simpler Hartree 
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approximation, with free space meson nucleon coupling constants = 8.79m^j_j;5 = 4.49m^ 



and g'^p = 3.86m^, with = 700 MeV, = 770 MeV and rup = 780MeV [ISj. Once we 
have the shell model core wave functions there is, however, no practical reason for not 
using the more sophisticated Hartree-Fock couplings for the hyperon. In a previous study 
of high central density neutron stars [l^, where the hyperon population is large enough 
that their exchange terms matter, we found that the Hartree-Fock couplings, g"^ = 11.33m^, 
g"^ = 7.27m^ and g^ = 4.56m^, gave a satisfactory phenomenology. So, for the hyperons we 
use these couplings (in Eqs.(4-7)). 

In the lsi/2 level of ^'^^Pb, this yields 26.9 MeV binding for a A and a mere 3.3 MeV 
binding for a The reduction of more than 23 MeV binding for the S° relative to the A is 
a tremendous improvement over the original QMC model and this may be traced directly to 
the inclusion of the effect of the medium in enhancing the hyperfine splitting of the S and 
A. It may be of interest to note that, if this hyperfine effect, which is an essential feature of 
the improved QMC model, were to be neglected, one would need, for example, to increase 
the u — T, coupling with respect to u; — A by about 20%. 

The remarkable agreement between the calculated and the experimental binding energy 
of the A in the lsi/2 level of ^osp^ jg another major success. In our earlier work the A 
was overbound by 12 MeV and we needed to add a phenomenological correction which we 
attributed to the Pauli effect. This correction is not needed when we use Hartree-Fock, 
rather than Hartree, coupling constants. 

TABLE II: Single-particle energies (in MeV) for y O, y Ca and y Ca hypernuclei. The experimental 



.a 



data are taken from Ref. [5l] (Table 11) for "O and from Ref. 



2M for 4"Ca. 





(Expt.) i^O 




^OCa (Expt.) 


A Ca 


4o Ca 




foCa 




-12.42 ±0.05 ±0.36 -16.2 


-5.3 


-18.7 ±1.1 


-20.6 


-5.5 


-21.9 


-9.4 


lP3/2 


-6.4 






-13.9 


-1.6 


-15.4 


-5.3 


lPl/2 


-1.85±0.06±0.36 -6.4 






-13.9 


-1.9 


-15.4 


-5.6 


11^5/2 








-5.5 




-7.4 




2si/2 








-1.0 




-3.1 




1^3/2 








-5.5 




-7.3 
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Already at this stage the binding of the E'' in the lsi/2 level of ^osp]^ jg j^g^- MeV 
- a major improvement over the earlier QMC results. However, as pointed out in that 
work, there is an additonal piece of physics which really should be included and which goes 
beyond the naive description of the intermediate range attraction in terms of a exchange. 
In particular, the energy released in the two-pion exchange process, NS^NA^NE, 
because of the E-A mass difference, reduces the intermediate range attraction felt by the 



E hyperon. In Ref. [25| this was modeled by introducing an additional vector repulsion for 
a E hyperon. Following the same procedure, we replace g^uj{r) by g^uj{r) + \y:Pb, with 
As = 50.3 MeV-fm^, as determined in Ref. 25f| by comparison with the more microscopic 



model of the Jiilich group 



32l |. We are comfortable using this earlier estimate of the effect 



of the coupled A — channel, even though there is a more recent potential from the Jiilich 
group 3J], because the latter tends to overbind hypernuclei 35| and the sign of the correction 



for coupling to a lower mass channel is in any case model independent. 

Our results are presented in Tables HTl and lllli The overall agreement with the experimen- 
tal energy levels of A hypernuclei across the periodic table is quite good. The discrepancies 
which remain may well be resolved by small effective hyperon-nucleon interactions which 
go beyond the simple, single-particle shell model. Once again, we stress the very small 
spin-orbit force experienced by the A, which is a natural property of the QMC model 26|. 

There are no entries for the E-hyperon because neither the E"^ nor the E° is bound to 
a finite nucleus. This absence of bound E-hypernuclei constitutes a major advance over 
earlier work. We stress that this is a direct consequence of the enhancement of the hyperfine 
interaction (that splits the masses of the E and A hyperons) by the mean scalar field in- 
medium. It is especially interesting to examine the effective non-relativistic potential felt by 
the E*^ in a finite nucleus. This is shown in Fig{T]for Calcium and Lead. In the central region 
the vector interaction dominates over the scalar one leading to a repulsive effective potential 
which reaches respectively 30 MeV and 12 MeV at the center. It is only at the surface that 
the scalar attraction becomes dominant. While the exact numerical values depend on the 
mass taken for the a meson, we stress the similarity to the phenomenological form found by 



Batty et al. 



361 ]. For a recent review see [37|]. It will clearly be very interesting to pursue 



the application of the current theoretical formulation to E~-atoms. 

We also note that this model supports the existence of a variety of bound S-hypernuclei. 
For the H° the binding of the Is level varies from 5 MeV in loO to 15 MeV in lo^Pb. The 
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TABLE II 
from Ref. 





fYh (Expt.) 


91 7„ 


91 Zr 


f^Ph (Expt.) 


^Pb 






-23.1 ±0.5 


-24.0 


-9.9 


-26.3 ±0.8 


-26.9 


-15.0 


lP3/2 




-19.4 


-7.0 




-24.0 


-12.6 


lPl/2 


-16.5 ±4.1 (Ip) 


-19.4 


-7.2 


-21.9 ±0.6 (Ip) 


-24.0 


-12.7 


1^^5/2 




-13.4 


-3.1 




-20.1 


-9.6 


2Si/2 




-9.1 








-17.1 


-8.2 


lc^3/2 


-9.1 ±1.3 (Id) 


-13.4 


-3.4 


-16.8 ±0.7 (Id) 


-20.1 


-9.8 


1/7/2 




-6.5 






-15.4 


-6.2 


2p3/2 




-1.7 








-11.4 


-4.2 


1/5/2 


-2.3 ±1.2 (1/) 


-6.4 





-11.7 ±0.6 (1/) 


-15.4 


-6.5 


2Pl/2 




-1.6 





— 


-11.4 


-4.3 


I99/2 




_ 


_ 


_ 


-10.1 


-2.3 


I97/2 








-6.6 ±0.6 (Ig) 


-10.1 


-2.7 


I/I1I/2 










-4.3 




2^5/2 










-5.3 




2^3/2 










-5.3 




I/I9/2 










-4.3 




3Sl/2 










-3.5 





experimental search for such states at facihties such as J-PARC and GSI-FAIR will be very 
important. 

IV. CONCLUDING REMARKS 

In conclusion, we stress that including the effect of the medium on the hyperfine inter- 
action between quarks within the quark-meson coupling model has led to some important 
advances. The agreement between the parameter free calculations and the experimental 
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FIG. 1: The effective potential for S*^ in Calcium and Lead 

ground state levels for A-hypernuclei from Calcium to Lead is impressive. In Oxygen the 
agreement is not as good and we take it as a reminder that the QMC model may be better 
at describing heavy nuclei, for which the constant density region dominates over the surface. 
For the d- and f-wave levels shown in Table III, there is a tendency for the model to overbind 
by several MeV. Whether this is a consequence of the use of an extreme single particle shell 
model for the core, the omission of residual A — interactions or an aspect of the current 
implementation of QMC that requires improvement remains to be seen. Nevertheless, we 
find these initial results, obtained with no adjustment of parameters to hyperon data, very 
encouraging. 

A number of S-hypernuclei are predicted to be bound, although not as deeply as in the 
A case. On the other hand, the additional repulsion arising from the enhancement of the 
hyperfine reulsion in the S-hyperon in-medium, together with the effect of the SA^ — AA^ 
channel coupling on the intermediate range scalar attraction, means that no S-hypernuclei 
are predicted to be bound. This encouraging picture of finite hypernuclei, suggests that the 
underlying model, which is fully relativistic and incorporates the quark substructure of the 
baryons, is ideally suited for application to the properties of dense matter and neutron stars. 
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